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O • Abstract. We present a coupled Boltzmann and hydrodynamics approach to 

QJ \ relativistic heavy ion reactions. This hybrid approach is based on the Ultra- 

1 relativistic Quantum Molecular Dynamics (UrQMD) transport approach with 

an intermediate hydrodynamical evolution for the hot and dense stage of the 
collision. This implementation allows to compare microscopic transport calcu- 
lations with hydrodynamic calculations using exactly the same initial condi- 

■ tions and freeze-out procedure. Here, the results of the different calculations 
qq , for particle multiplicities, the mean transverse mass excitation function and 

the HBT radii at all SPS energies are discussed in the context of the available 
. data. 
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One of the main motivations to study high energy heavy ion collisions is the 
' creation of a new deconfined phase of strongly interacting matter, the so called 

Quark-Gluon Plasma (QGP) [[lj. Since the direct detection of free quarks and 
gluons is impossible due to the confining nature of QCD, it is important to model 
the dynamical evolution of heavy ion reactions to draw conclusions from the final 
state particle distributions about the interesting early stage of the reaction. To get 
a more consistent picture of the whole dynamics of heavy ion reactions various so 
called microscopic plus macroscopic (micro+macro) hybrid approaches have been 
launched during the last decade [HJ [3 [H [H [6l H [HI [^). 

In this paper we briefly describe the specific micro+macro hybrid approach that 
embeds a hydrodynamic phase in the UrQMD approach. This allows to reduce the 
parameters for the initial conditions and the freeze-out prescription. At present cal- 
culations imposing a hadron gas EoS are shown to provide a baseline to disentangle 
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the effects of the different assumptions for the underlying dynamics in a transport 
vs. hydrodynamic calculation. We show results for particle multiplicity, HBT radii 
and (tot) excitation functions in the context of the available data at FAIR/SPS 
energies. 

The Ultra-relativistic Quantum Molecular Dynamics Model (UrQMD) is used 
to calculate the initial state of a heavy ion collision for the hydrodynamical evolution 
[ O QUI [TT] . This has been done to account for the non-equilibrium dynamics in 
the very early stage of the collision. In this configuration the effect of event-by- 
event fluctuations of the initial state are naturally included. The coupling between 
the UrQMD initial state and the hydrodynamical evolution proceeds when the two 
Lorentz-contracted nuclei have passed through each other. 
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Fig. 1. (Color online) Excitation function of particle multiplicities (4tt) in 
Au+Au/Pb+Pb collisions from Ei ab = 2 A GeV to y/IWW = 200 GeV. 
UrQMD+Hydro (HG) calculations are depicted with full lines, while UrQMD-2.3 
calculations are depicted with dotted lines. The corresponding data from different 
experiments [[ISl[I71[Il[Il[ini[IIl[ll[ll[ll[lSl[ia[17]are depicted with 
symbols. 
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After the UrQMD initial stage, a full (3+1) dimensional hydrodynamic evolu- 
tion is performed using the SHASTA algorithm [[TJ1 [T3]. For the results presented 
here an equation of state for a free hadron gas without any phase transition is used 

The hydrodynamic evolution is stopped, if the energy density drops below five 
times the ground state energy density (i.e. ~ 730MeV/fm 3 ) in all cells. This cri- 
terium corresponds to a T-//B-configuration where the phase transition is expected 
- approximately T = 170 MeV at \ib = 0. The hydrodynamic fields are mapped 
to particle degrees of freedom via the Cooper-Frye equation on an isochronous hy- 
persurface. The particle vector information is then transferred back to the UrQMD 
model, where rescatterings and the final decays are perfomed using the hadronic 
cascade. A more detailed description of the hybrid model including parameter tests 
and results for multiplicities and spectra can be found in [ 1 1 5 J . 

In the following results from the UrQMD+Hydro (HG) model are compared to 
the default transport calculation (UrQMD-2.3) and to the experimental data. Note 
that parameters are constant for different energies and/or centralities in the hybrid 
model calculation. _ , 
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Fig. 2. (Color online) Mean transverse mass excitation function of pions at 
midrapidity (|y| < 0.5) for central (b < 3.4 fm) Au+Au/Pb+Pb collisions from 
Eub = 2 - 160 A GcV. UrQMD+Hydro (HG) calculations are depicted with full 
lines, while UrQMD-2.3 calculations are depicted with dotted lines. The corre- 
sponding data from different experiments [[28l HH [20] are depicted with symbols. 

Fig. Q] shows the excitation functions of the total multiplicities for central 
Au+Au/Pb+Pb collisions for E\ a \> = 2A GeV to ^/snn = 200 GeV. Compared 
to the default simulation, the pion and proton multiplicities are decreased over 
the whole energy range in the hybrid model calculation due to the conservation 
of entropy in the ideal hydrodynamic evolution. The non-equilibrium transport 
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calculation produces entropy and therefore the yields of pions are higher. The 
production of strange particles however, is enhanced due to the establishment of 
full local equilibrium in the hybrid calculation. Since the yield of strange particles 
is small they survive the interactions in the UrQMD evolution that follows the 
hydrodynamic freeze-out almost without re-thcrmalization. 

Next, we turn to the excitation function for the mean transverse mass of pions. 
The UrQMD approach shows a softening of the equation of state in the region where 
the phase transition is expected because of excited resonances and non-equilibrium 
dynamics while the hydrodynamic calculation with hadron gas EoS just rises as a 
function of the energy. 
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Fig. 3. (Color Online)fcx dependence of the Pratt-radii of the pion source in central 
Pb+Pb collisions at £q ab = 20, 30, 40, 80, and 158A GeV and mid-rapidity < 
Y-k-k < 0.5. Preliminary NA49 data are taken from [I30j. 



Finally we show HBT results. It is well known that the Hanbury-Brown- 
Twiss interferometry (HBT) technique can provide important information about 
the spatio-temporal structure of the particle emission source (the region of ho- 
mogeneity) [ [29]. The HBT radii can be extracted by fitting the correlator of 
the particle pair with a Gaussian form. Using the Pratt's three-dimensional con- 
vention (the LCMS system), the parametrization of the correlation function reads 
C(qL,qo,Qs) = 1 + Xe^ R2Lq2L ~ R oio- R s'Js- 2R oLiogL w here A is normally referred to 
as an incoherence factor, but we can regard it as a free parameter since it might be 
affected by many other factors. Rl, Ro, and Rs are the Pratt radii in longitudinal, 
outward, and sideward directions, while the cross-term Rql plays a role at large 
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rapidities, qi is the pair relative momentum q = Pi — P2 in the i direction. 

Fig. [3] shows the transverse momentum fcx dependence (kj 1 = (Pit + P2t)/2) 
of the HBT- radii Rl (top plots), Ro (middle plots), and R$ (bottom plots) of the 
negatively charged pion source in Pb+Pb reactions at beam energies £i a b = 20, 30, 
40, 80, and 158A GeV. The calculations (lines with symbols) are compared with 
preliminary NA49 data (stars) [ 30J (< 7.2% of the total cross section <tt)- The 
pair-rapidity < < 0.5 is chosen for all reactions. 

One observes that the calculations in the UrQMD cascade mode can reproduce 
the fey-dependence of HBT radii Rl and Rs fairly well. Only at small fcx values, 
the calculated Rl and Rs values are up to 25% lower than data [[31]. For the Ro 
values, the calculations reach higher values than the experimental data especially at 
relatively large fey. The HBT Rl radii from the hybrid model are seen to be larger 
than data and those from cascade UrQMD calculations. Interestingly, however is 
the fact that the stronger expansion leads to a better description of the R s radii, 
improving the R Q / R s ratios. Certainly, the EoS used in hydro-process in the hybrid 
model will also strongly influence the final HBT radii and will be investigated in a 
forthcoming paper. 
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Fig. 4. (Color online) Beam-energy dependence of the ratio Ro/Rs at 100 < kr < 
200 MeV/c for central Pb+Pb collisions at SPS energies. 



Fig. d] illustrates the beam-energy dependence of the ratio Ro /Rs of the pion 
source at 100 < kx < 200 MeV/c for central Pb+Pb collisions at SPS energies. In 
the cascade mode, the calculated ratio is higher than the experimental data, while 
with the embedded hydrodynamic evolution the ratio is driven down towards the 
experimantal data. With increasing beam energy, the difference of the Ro / Rs ratio 
between the default transport and the hybrid model calculation becomes larger. 

We have presented the main ideas of an integrated Boltzmann and hydrody- 
namics approach to relativistic heavy ion reactions. The final pion and proton 



G 



H. Petersen et al. 



multiplicities are lower in the hybrid model calculation due to the isentropic hy- 
drodynamic expansion while the yields for strange particles are enhanced due to 
the local equilibrium in the hydrodynamic evolution. The results of the different 
calculations for the mean transverse mass excitation function and the HBT radii 
at all FAIR/SPS energies are sensitive to the underlying dynamics and reflect the 
stronger transverse expansion in the hybrid model. 
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